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Abstract: The structural, electronic, optical and thermoelectric characteristics of crystalline 
oxides-perovskites BaTMOs were investigated using all-electrons full-potential linearized 
augmented plane wave (FP-LAPW) method within the framework of density functional theory 
(DFT). These calculations are performed using the computer code WIEN2k. The generalized 
gradient approximation as parameterized in Perdew, Burke and Ernzerhof (PBE-GGA) is 
employed to calculate exchange-correlation potential. Also, the modified Becke Johnson 
exchange potential approximation as parameterized by Tarn and Blaha (TB-mBJ) is used to 
improve the band gap estimation. The comparison among the previous ab intio calculations, the 
experimental measurements, and our study results are done. In our calculations, it is shown that 
these two perovskites BaZrO3 and BaHfOs are insulators and have widely indirect band-gap 
energy (К-Г) from both approaches, РВЕ-ССА and TB-mBJ. The optical properties such as 
dielectric tensor, the refractive index, the absorption coefficient, and the electron loss function 
have been calculated and analyzed. Furthermore, the transport characteristics calculations based 
on semi-classical Boltzmann theory have been discussed. 


Keywords: Perovskite, DFT, PBE-GGA, TB-mBJ, Optical, Ferroelectric, BaZrO3, BaHfOs, 
Thermoelectric. 


*corresponding author (Said Al Azar; Email: q_saed74@yahoo.com) 


1. Introduction 


Both Barium zirconate (BaZrOs) and Barium hafinate (BaHfO3) perovskites get more 
attention experimentally ‘and theoretically 8 '3 because of their unique physical properties and 
potential applications as functional materials. They show an extensive variety and fascinating 
such as high thermoelectric power, ferroelectricity, piezoelectricity frequency doublers, 
superconductivity, giant magnetoresistance, and the exchange of structural, optical, magnetic, 
electronic characteristics 21515, The Perovskites (ABX3), A and B are cations and X is an anion, 
there are many kinds of the perovskite structure crystals, such as the oxides-perovskites (ABO), 
flouro-perovskites (ABF3), anti-perovskite (A3BX), double-perovskite (A2BOe) and nitrides-anti- 
perovskites (A3BN)!6. 


Barium zirconate, BaZrOs is a surface effect ceramic material that has many mechanical 
and technological applications". Also, BaZrOs has some unique properties such as high 
dielectric constant (-50). And it has high thermal stability and wide bandgap energy (-4 eV). 
These diverse properties have placed this perovskite ceramic as a useful material in applications 
like microwave devices, and high-temperature proton conductor for electric ceramics ?. King- 
Smith and Vanderbilt calculated lattice constants, band structures, zone-center phonon 
frequencies, and elastic constants using the FLAPW method, with the exchange local density 
approximation (LDA) 5. Barium hafnate BaHfOs, is the compound have high-temperature 
melting (2893K), physical unique properties have been studied in both experimental and 
theoretical terms including thermal expansion coefficient, heat capacity, thermal conductivity, 
hardness, by using the X-ray diffraction pattern, and the lattice parameter of 0.4171 nm at room 
temperature 18, Also, Vali studied the properties such as the electronic band structure, phonon 
dispersion and dielectric constants of BaHfOs using DFT. He found that BaHfOs has an indirect 
bandgap between (В-Г) points about 3.533 eV”. 


Thin films of ВАНЬ х ТІОз are deposited using pulsed liquid injection MOCVD technique on 
Si (100) substrate. The roughness average value (Ra) of the film is increasing with increase x- 
value. Also the BaHfo.7sTio.25O3 thin film shows well-behavior capacitance-voltage characteristic 
and its k-value is around 14 20. Evarestov investigated the structure stability for both BaZrOs and 
BaHfOs perovskites using hybrid functional DFT with phonon frequency calculations. The 
Mulliken and Born atomic charges results are emphasized the stability of cubic and covelant 
bonds іп both crystals 12. 


Photoluminescence and photocatalytic activity of bismuth doped BaZrOs were investigated 
experimentally using powder X-ray diffraction (XRD), UV-VIS diffuse reflectance 
spectroscopy, photoluminescence (PL) spectroscopy, and photocatalytic activity. The results 
suggested that BaZrOs:Bi could be a good candidate to be used as a visible-light activated 
photocatalyst under excitation wavelengths less than 800 nm ?!. Parida et al synthesized BaZrO3 
ceramic powder by the solid state reaction method. Their measurement indicated that BaZrOs 
could be used for small dielectric resonator antenna (DRA) technology because of its low less 
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microwave dielectrics materials ?. Khan and Qureshi observed that Tantalum doped BaZrOs at 
Ta concentration х = 0.04 (BaZro.96Tao.0403) gives the highest rate of hydrogen evolution 
(180umol/h) along with the surface area 23. 


Zhang et al conclude that the spontaneous polarization, Born effective charge and dielectric of 
BaZrOs is significantly sensitive to epitaxial strain and the transition from paraelectric to 
ferroelectric happened in both compressive and tensile strain. BaZrOs could be useful for tunable 
dielectric materials >. Chen et al synthesis BaZrOs hollow nanostructure with oxygen vacancies 
through one step solvothermal method using ethylenediamine (EDA) as solvent for 
photocatalytic hydrogen evolution from pure water. The H» production rate is about 15 times 
than that of the commercial BaZrO3 (BZO-C) ?. The BaZrO3 hollow nanocrystal dopant by Бе"? 
ions were prepared through facile solvothermal method for photocatalytic Н> from pure water 
under visible light. The existence of Fe (Ш) ion plays a key role in broadening the optical 
absorption range and restricting the recombination of photo-induced carriers under visible light 
irradiation 852027, 


Although these two perovskites are studied intensively, the literature review indicates that 
there is a paucity of thermoelectric data and calculations. The main goal of this research study is 
to review and investigate in details the ferroelectric and thermoelectric properties of the cubic 
Barium Zirconate and Barium Hafnate perovskites to design new potential applications in 
ferroelectric and thermoelectric technologies. 


This contribution is organized as follows: an introduction with the previous overview, the 
computational and theoretical framework details, the results and discussion with the comparison 
with previously available data, and finally the conclusions are drawn. 


2. Computational Details 


We have performed total energy calculations for the oxides-perovskite ВаТМО», an ideal 
cubic structure that contains only one, composed of four atoms in its unit cell. The cations Ba are 
at the cube corners, Та (0, 0, 0), while the cations TM are at the center of the cube, 1b (1, 12, 12) 
and oxygen atoms at the center of the cube faces forming a regular octahedron, Зс (0, >, 12) sites 
of Wyckoff coordinates, as shown in Figurel. The space group for these two perovskites is Pm- 
3m (#221). We have used the interstitial region to express plane waves are used approximation 
for the calculation of exchange-correlation energy functional, namely; generalized gradient 
approximation (GGA) of Perdew, Burke and Ernzerhof (PBE) parameterization 28 and using 
modified Becke Johnson exchange potential approximation as parameterized by Tran and Blaha 
(TB-mBJ) to calculate the electronic band structure and density of states 99. The TB-mBJ 
method was known that it provides more accurate energy band gap values than PBE-GGA 
method with cheap calculations cost. Where the first-principles calculations were performed 
based on (DFT) to describe the ground state ? used the method of (FP-LAPW) 2355 this is 


performed using the computer code іп Wien2k code 39. Also, we chose the interstitial area of the 
plane wave with a cutoff Клах X Кит = 8, where Рут is the smallest muffin-tin radius in the 
unit cell, also Ктах is the largest К vector magnitude of the interstitial plane-wave expansion of 
the basis set. The chosen number of k-points is (16x16x16) in codes, and which corresponds to 
the total number of k-points of 4096. 


For the optical properties, we calculated the complex dielectric function ¢(@) = єго) + 1€2(@), 
where &1(0) is the real part and £»(o) is the imaginary part. The real part ¢1(@) of dielectric 
function e(w) gives from the Kramer-Kronig relationship: 





2 оо WE2(0) ; 
в1(0) = 1+ = Р A ура do [1] 
where the P is the principle value of the integral. The imaginary part ¢2(@) was calculated from 
the momentum matrix elements between the occupied and unoccupied wave functions within the 
selection rules 
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where р is Ше momentum matrix element between states of bands a and В with crystal 
momentum К. ск and ук are the crystal wave functions corresponding to the conduction, and 
valance bands with crystal wave vector k 37. There are two contributions to e(w), namely; inter- 
band and intra-band transitions. The contribution inter-band is important only for metals, where 
be inter-band between valence-to-conduction absorption spectra. Whereas intra-band between 
valence-to-valence and conduction-to-conduction absorption. Which is describes the complete 
response of a material to the applied electromagnetic radiations field 28. Then, the main optical 
parameters: the refractive index n'‘(w), the extinction coefficient k“(w), the absorption 
coefficient a" (а), optical conductivity oa" (w), electron loss energy function LË (w), and the 


normal incidence reflectivity RË (w) in the crystal are given by 240. 


nii Gre јез (а) + £l? (0) + lw) 
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Е (о) = [8] 
The transport properties are obtained from Ше analysis of band structure results based on 
Boltzmann theory method under the relaxation time approximation as implemented in BoltzTraP 
code 72, To explore the properties of thermoelectric such as Seebeck coefficient (S), electrical 
conductivity (o), specific heat (Ce), electronic thermal conductivity (x^), and power factor rate 
(PF). A dense k-mesh of 25000 in the full Brillouin zone was used to guarantee good results. 


The Seebeck coefficient (S) is related to carrier concentration via the Mott formula as follows *: 


_ КВТ 1 ап(Е) | таи(Е) 
S= 3e 5 dE и ти [9] 


where е is the electron charge, Кв is Boltzmann constant, h is the Planck constant, u is carrier 
mobility and п is the carrier concentration. 


Electrical conductivity (с) is related to carrier concentration as follows": 


пет 


о = [10] 


т 





where m is Ше electron mass, п is the carrier concentration апа т is Ше relaxation time. 


Electronic specific heat (Ce) is related to temperature and chemical potential as follows +: 


d Т; 
Ce(T; н) = fnt) - и) [| ae [11] 
3. Results and Discussion 


3.1 Structural Properties 
Goldschmidt tolerance factor (t) for perovskite is given Бу“: 


— (ratrx) 
С V2(rg+rx) [12] 


Where ryis the ionic radius of atom N. The t values of existence perovskite compounds range is 
(0.75 « t € 1). For ideal cubic perovskite structure, t has unity value. The tolerance factor t for 
BaZrO; and BaHfOs are 1.0 and 1.016, respectively. 


We presented the relationship between the energy changes as a function of the unit cell 
volume. Using the Birch-Murnaghan equation of state to determine the different important 
structural and ground-state properties “7. The equilibrium lattice constants (а) and ground-state 
energies (Eo) are calculated, as estimated from optimization at zero pressure. Besides the volume, 
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Ше lattice constants increase because of Ше increase in Ше atomic radius from Zr to Hf. In Ше 
present work, we were fitting the energy-volume values to a third-order Birch Murnaghan 
equation of state to calculate lattice parameter a and bulk modulus Bo of Ba(Zr, Hf)O3 
compounds, where the bulk modulus increases as the volume decreasing with moving of cation 
down in the group in the periodic table. The calculated ground-state parameters for both BaZrOs 
and BaHfOs are tabulated in Table 1 along with previous computational and experimental 
studies. For both perovskites, our results are in reasonable agreement with the available 
experimental data and previous calculations. 


3.2 Electronic properties 


То study the electronic properties of BaZrO3 and BaHfO3, DOS’s and band structures along 
the symmetry direction of the first Brillouin zone, of these compounds both PBE-GGA, and TB- 
mBJ methods for exchange-correlation potential are used. The band structures for BaZrOs and 
BaHfOs are presented in figures 2 and 3 respectively. The Fermi level Er is chosen to locate at 0 
eV. In both methods, the BaZrOs shows semiconductor characteristics with indirect bandgap (Eg) 
(В > Г), but in the TB-mBJ method, the (Ez) is wider. The valence band maximum (VBM) 
occurs along the R-point symmetry line, while the conduction band minimum (CBM) occurs 
along the Г-роіпі symmetry line. The bandgap approaches 3.39 eV by using the PBE-GGA 
method whereas it approaches 4.424 eV by using the TB-mBJ method. As well, BaHfO3 shows 
the same behavior with indirect semiconductor bandgap (R -» Г), where the bandgap 15 3.69 eV 
by using the PBE-GGA method and 5.246 eV when using the TB-mBJ method. Table 2, the 
calculated bandgap (Eg) of both perovskites BaZrOs апа BaHfO3 are compared with previous 
theoretical and experimental data. For both compounds, the TB-mBJ values approach more 
enough to experimental results. 


The total density of state (TDOS) and partial density of state (PDOS) site projected for both 
compounds are calculated using both methods PBE-GGA and TB-mBJ. The TDOS and PDOS 
for BaZrOs are shown in figure 4. It is clear that the lower valence band consist of Ba-p state 
around -10.1 to -11 eV, and the upper valence band extends down to approximately -4 eV and is 
composed mainly of O-p states with a small contribution from Zr-d states, while the lower 
conduction band is mainly a mixture of Ba-d and Zr-d states, with a small admixture from O-p. 
In the conduction bands, the Zr-d states are dominant. Figure 5 displays the TDOS and PDOS for 
BaHfO3 perovskite using both PBE-GGA and TB-mBJ methods. The contribution to the lowest 
valence band at - 10 eV 15 due to Ba-p state, whereas the upper valence band (around -4.2 to 0 
eV) is from O-p states hybridized with some Hf-p and Hf-d electrons. The bands above the 
Fermi level are the conduction band of the compound; the lower part of the band (around 3.8 to 6 
eV) is dominated by the O-p. Also, it can be seen from the plot, the top of consist are mainly Ba- 
d and Hf-d states hybridized with some O-p. The conduction bands are composed mostly of 
antibonding between Hf-d and O-p as well as antibonding between Ba-d and O-p, maximum 
peak at 7.5 eV 15 due to the unoccupied d states of Hf. 


3.3 Optical Properties 


In this subsection, the optical characteristics of BaHfOs and BaZrOs are presented. 
Electronic excitations will take place within the crystal if it is exposed to electromagnetic 
radiation. As a result, the optical features of the crystal will be observed. The nature of the 
response to an incident electromagnetic wave at a particular frequency determines if this material 
could be used in optoelectronic devices. The different optical characteristics of the crystals are 
mainly calculated by using the frequency-dependent complex dielectric tensor #(а)). Since 
BaZrOs and ВаНіОз perovskites has cubic structure symmetry which means one tensor 
component is enough to do optical calculations. Thus, the optical characteristics of the materials 
n(o) , k(o), R(w) , o (o) , a(w), and L(w) could be derived from their =; (w) and €2(w) values 
as illustrated in the computational details section +. 


The imaginary part €,(w) is determined from the electronic structure through the joint 
density of states and the momentum matrix elements between the occupied and the unoccupied 
wave functions within the selection rules. In figure 6, the optical characteristics for both BaZrOs 
апа BaHfOs compounds are plotted. In panel a) the spectrum of €,(w) shows the three peaks 
structure for BaZrO3 and BaHfO3. The є» (о) has deeply impact on the absorption spectrum of 
the matter. A large value of £5(«) refers to a large absorption spectrum include of the whole of 
all different transitions that will occur between the valance and the conduction bands. It is 
starting to response to the incident radiation at critical points which are the threshold energy at 
2.87 eV апа 4.26 eV for BaZrOs and BaHfOs respectively. These points are Cv—Cc splitting and 
give a threshold for indirect optical transitions between the highest valance band and the lowest 
conduction band. This is known as the fundamental absorption edge. The main peak in the 
spectrum is situated at 7.12 eV and 7.52 eV for BaZrOs, and BaHfOs, respectively. These peaks 
are dominated by transitions from the O-p band and the small contribution of transition metals d 
state just below Fermi energy to Ba-p state of the conduction band. In panel b) the spectrum of 
є (00) for BaZrO3 and BaHfO3 compounds show that the static dielectric constant ғ1(0) are 
equal 3.69 and 3.51, whereas the main peaks appear at energy 5.24 eV and 5.59 eV, respectively. 
It is important to consider both sı and e2 together since they affect each other, meaning the shape 
of є» cause corresponding changes in the shape of є; and vice-versa. This is known as the 
Kramers-Kronig relation between the real (£1) and imaginary (22) parts of the dielectric function. 
In short, the dielectric function describes what an electric field such as an oscillating light wave 
does to the material. In panel c) the absorption coefficients for BaZrO3 and BaHfO3 compounds 
show the same behavior as є. (w). The absorption edges start from 3.28 and 4.29 eV for BaZrOs 
and BaHfOs respectively which are corresponding to optical bandgaps. This produces from the 
transition of the O-p valence band electrons to the empty d conduction bands in Zr or Hf. In 
ВаНГО: perovskite, the main peak in dielectric absorption spectra occurs at 7.81 eV whereas the 
two other peaks at 11.3 and 11.78 eV whereas in BaZrOs the main peak occurs at 7.3] eV the 
other two peaks occur at 8.8 and 11.0 eV. In both compounds, the Ba does not contribute to the 
low-energy absorption spectra; because its electron valence bands are not in the transition energy 


domain. In panel d) Ше loss function for BaZrOs and BaHfO3 compounds are plotted and their 
peaks occur at 7.6 eV and 8.4 eV, respectively. In panel e) the reflectivity for both BaZrO3 and 
BaHfO3 compounds are plotted and their static reflectivity R(0) is approximately equal 10%. At 
low energies (1.е., 0—7.8 eV for BaHfO3 and 0-7.2 eV for Ва/гОз), the reflectivity spectrum was 
observed to extend as high as 33% and 30% for BaHfO3 and BaZrOs;, respectively. Our 
calculations show a strong reflectivity at high energies increase up to 75%, and 72% for BaHfOs, 
and BaZrOs, respectively. In panel f) the refractive index plots show that the static refractive 
index n(0) are equal 1.92 and 1.87 for BaZrOs and BaHfO;, respectively. The refractive index 
reaches its maximum value 2.7 at 7.5 eV for BaHfOs, whereas for BaZrO3 its 2.5 occurs at 5.0 
eV. 


Finally, the optical conductivity o(w) for both BaZrOs and BaHfOs perovskites are displayed 
in figure 7. There is quite similarity in curves of both perovskites. For BaZrOs the first peak 
appears at 6.8 eV, the second peak appears at 9.0 eV whereas the third peak at 10.5 eV. On the 
other hand, at BaHfO3 there are two mean peaks at 7.4 eV and 11.5 eV. The maximum value for 
BaZrOs is equal 8900 at 10.5 eV while for BaHfOs it is equal 9100 at 11.5 eV. 


3.4 Thermoelectric (TE) Properties 


The efficient thermoelectric characteristics for any materials determines by their figure of 
merit zT = 520Т/к°, the increasing of S and c increase figure of merit while increasing к? 
decrease figure of merit. As we known, a suitable material for thermoelectric applications should 
have figure of merit equal or greater than one which give energy transformation efficiency 
greater than 25% ^9. 

Figure 8 shows the variation of Seebeck coefficient S with respect to temperature for both 
compounds BaZrOs and BaHfO3. The two curves have the approximately same behavior where 
they change nonlinearly by increasing from zero to reach the maximum value 206.86 иУ/К at 
T=270K for BaZrO3 and 288.13 uV/K at T=180K for BaHfOs, after that it is decreasing slightly 
in BaHfO3 but in BaZrOs it is still increasing. At room temperature (RT), Seebeck coefficients 
for BaZrO3 and ВаНЮз are 208 and 261 цУ/К, respectively. In figure 9, the electric 
conductivity c as a function of temperature when the constant relaxation time то = 0.8 x 1075 
is used, it is shown that the electric conductivity is decreasing non-linearly at the low 
temperature to reach its global minimum at T=121K in BaZrOs case and Т=90 in BaHfOs, then it 
starts increasing linearly with temperature increases for both BaZrO3 and ВаН Юз. The electric 
conductivity at RT for BaZrOs is about 0.554 x 1055 / т and it about 0.289 x 1055 / т for 
BaHfO3. Figure 10 displays the electronic thermal conductivity x? as a function of temperature 
when the constant relaxation time ту = 0.8 x 107175 is used. The behavior of the electronic 
thermal conductivity curves are nonlinear increasing for both compounds. The electronic thermal 
conductivity of Ван Оз at RT is equal to 0.154 Ма КО, which is about half the BaZrO3 value 
of 0.309 Му КЛ. Our results contradict the previous experimental values 10.4 ?? and 5.22! Wnr 
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IK", respectively. This contradicts maybe because we calculated only the electronic part of 
thermal conductivity but the experiment results were for total thermal conductivity where the 
lattice part was dominant in both compounds. The temperature dependence of the thermal 
conductivities of Ва/гОз and BaHfO3 follow an inverse temperature law (1/T) at low 
temperatures, showing typical phonon conduction characteristics ?!. 


Figure 11 shows how are the power factor (PF) curves changed when the temperature 
increase for both compounds. They reach their maximum values 0.213 un W/mk? at T=180K for 
BaHfOs and 0.0814 uW/mk? at Т-240К for BaZrOs. Also, at RT the PF's equal 0.128 uW/mk? 
for ВаН Оз and 0.0778 uW/mk? for BaZrOs. Although the PF values for both compounds аге 
very small if they compared with the well-known thermoelectric materials but the results show 
that BaHfOs is better than BaHfOs if a suitable dopant used. 


To determine in details thermoelectric (TE) behavior for both BaHfOsand BaZrO;, such as 
Seebeck coefficient (S), electrical conductivity (с), power factor (PF), electronic specific heat 
(Се), and electronic thermal conductivity (к') have been expressed as a function of chemical 
energy (и); the calculated characteristics were plotted for considered temperatures: 300 - 800 K. 
In figure 12, the Seebeck coefficients S for both BaZrOs and BaHfOs compounds as a function 
of chemical energy и at different temperature are plotted. For both compounds, the curve 
increases linearly from zero value to reach the maximum value and then decreases sharply to 
negative global minimum values after that it increases linearly again to reach zero value. 


The chemical potential dependencies of electric conductivity rate (с/т) for both BaZrOs and 
ВаНЮз compounds at different temperatures are depicted in figure 13 panels а) and b). In the 
two panels, the electric conductivity goes to zero in energy range 0.52 — 0.74 Ry for BaZrOs and 
0.54 — 0.8 Ry for BaHfOs for the wide range of temperatures from 200K to 800K. In figure 14 
the power factor with respect to chemical energy at different temperature are plotted for BaZrOs 
and ВаНЮз perovskites. The two plots display two main peaks 11.7 x10!! W/mK?s and 6.0 
x10!! W/mK?’s at 0.525 Ry and 0.83 Ry for BaZrOs respectively, and 10.65 x10!! W/mK?s and 
6.75 x10!! W/mK?s at 0.5 Ry and 0.775 Ry for ВаН Оз respectively. Figure 15 shows the 
electronic specific heat capacity as a function of chemical energy at different temperature. 
Cordfunke and Konings have estimated the heat capacity at room temperature from a comparison 
with BaZrOsand recommended 101.7 ТК тої ‘at 298.15 К >. The calculated heat capacity of 
BaZrOs in this study is well agreement with the data in the literature 2273, Figure 16 shows the 
variation of the density of state around the Fermi level with respect of temperature. For the two 
alloys BaZrOs and BaHfOs the band gap still appear for a wide range of temperature. 


4. Conclusions 


Employing the all-electrons full-potential LAPW method based on DFT, with using PBE-GGA 
functional of exchange-correlation energy; and semi-local exchange potential TB-mBJ to study 


physical properties. Also, semi-classical Boltzmann transport theory has been used. The 
structural, electronic, optical properties, as well as the transport properties of BaZrOs and 
BaHfOs;, have been studied. A summary of results are: 


(1) The calculated lattice constants at zero pressure of these compounds are in a reasonable 
agreement with the previous theortical and experimental values. 


(ii) The electronic bandstructure calculations showed that BaZrO3 and BaHfOs are insulators 


materials with an indirect bandgap (А – Г). The fundamental bandgap of this compounds 
increases when mBJ-GGA approach is used instead of PBE-GGA. 


(iii) The joint DOS and the Kramers-kronig transformation are used to calculate and analyze the 
optical characteristics: dielectric constante(w), the absorption coefficienta(w), the 
reflectivityR(w), the refractive index n(w) and electron loss energyL(w). The oxygen p states 
and Transition metals Zr- and Hf-d states play a major role in optical transitions. 


(iv) Finally, The calculated transport characteristics showed that the two compounds are not 
suitable choice for thermoelectric applications but if they dopant with some impurities their 
properties could be improved. 
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ТаМе1: Lattice constant a(A), bulk modulus B(GPa), for the perovskties BaZrO3 and BaHfOs. 








Compound a(A) B(GPa) Method 
4.234 149.8 This work 
4.1895, 4.1935, 4.196855, 4 19245156. 1035156 155.910, Exp 
S 4.226910, 4.206377, 4.181222, 4.1815, 127” 
4.19327 
4.1545 4.207", 4.14818, 4.2560, 4.19313, 175.38,157, 174.718, | Theory 
4.162851, 4.192, 4.19812 1415, 1722, 1682 
4.212 156.0 This work 
4.171520, 4.1762, 4.17263, 4.16686, 4.179655, Ехр 
зн 4.17113, 4.192, 4.19312, 4.31066, 4.24991, 14413, 1752, 1712, Theory 
4.170P, 4.286", 4.2489", 4.16088, 4.0649? 123.695. 1869. 
19411, 146.721, 
145.50 ^*. 17969 








Па, Ма Obtained from the Voigt-Reuss-Hill (VRH) approximation. 


11b, 14 Obtained from fitting the third-order Birch-Murnaghan equation 
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Table 2: The calculated bandgap (Eg) using PBE-GGA апа TB-mBJ methods with comparison 
with previous computational and experimental results 








Ва2гОз BaHfO; 
Eg (eV) method Ref. Е, (eV) method Ref. 
4.42 TB-mBJ This work 5.25 TB-mBJ This work 
3.39 PBE-GGA This work 3.69 PBE-GGA This work 
5.3 Exp. р 6 Ехр P 
3.9 LDA 18 2.99 LDA 11 
3.156 ССА v 3.94 SX-LDA It 
5.4 LCAO 32 3.17 GGA is 
3.55 Exp. 2 5.3 SX-LDA 69 
4.96 Ехр. 28 3.9 ССА-РВЕ Pd 
5.9 TB-mBJ 58 





Figure Captions: 


Figurel: Crystal structures of BaTMOs (TM: blue, Ba: green and O: red). 

Figure2: Band structure of cubic perovskite BaZrOs using a) PBE-GGA and b) TB-mBJ 
methods. 

Figure3: Band structure of cubic perovskite BaHfOs using a) PBE-GGA апа b) TB-mBJ 
methods. 

Figure4: Calculated total and Partial-projected DOS for Ва/гОз using PBE-GGA compared 
with TB-mBJ method. 

Figure5: Calculated total and Partial-projected DOS for BaHfO3 using PBE-GGA compared 
with TB-mBJ method. 

Figure6: Calculated optical properties for the BaHfO3 and BaZrOs compounds, а) ¢2(@), b) 
£1(0), с) a(@), d) (о), e) R(o), and f) n(o). 

Figure7: Calculated optical conductivity o(@) for BaHfO3 and BaZrO3 

Figure8: Seebeck coefficient as a function of temperature at ground state chemical potential for 
the BaHfO3 and BaZrO3 compounds. 

Figure9: The electrical conductivity rate (с/т) as a function of temperature at ground state 
chemical potential for the BaHfO3 and BaZrO3 compounds. 

Figure10: The electronic thermal conductivity x? as a function of temperature at ground state 
chemical potential. 

Figurell: The power factor (PF) as a function of temperature. 

Figure 12: Seebeck coefficient S as a function of chemical potential at different temperature for 
a) BaHfO3 and b) BaZrOs. 

Figure13: (a) and (b) are the electrical conductivity rate (с/т) as a function of chemical potential 
at different temperature for both BaHfO3 and Ва/гОЗз; (c) and (d) are the electronic thermal 
conductivity as a function of chemical potential at different temperature for both compounds 
BaHfOs and BaZrOs, respectively. 

Figureld: The power factor (PF = 52с/т) as a function of chemical potential at different 
temperature for а) BaHfO3 and b) BaZrOs. 

Figurel5: The electronic specific heat as a function of chemical potential at different 
temperature for a) BaHfOs and b) BaZrOs. 

Figure 16: The density of state at Fermi level DOS[Er] as a function ої chemical potential at 
different temperature for a) BaHfO3 and b) BaZrOs. 
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Highlights (for review) 


Highlights: 


1- In this study, the electronic, optical and thermoelectric properties of Ba(Zr, Hf)O3 
perovskites are investigated computationally using LAPW method as implemented in 
DFT. 


2- The two compounds exhibit insulator properties with wide indirect band gap. 


3- The two compounds have good optical properties that make them potential candidates 
for optoelectric applications. 


4- The transport properties for both compounds are not as good to use in thermoelectric 
applications. 
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1. Introduction 


Both Barium Zirconate (BaZrO3) and Barium Hafnate (BaHfO3) perovskites get more 
attention experimentally" and theoretically" !! because of their unique physical properties and 
potential applications as functional materials. They show an extensive variety and fascinating 
characteristics such as ferroelectricity, piezoelectricity frequency doublers, superconductivity, 
giant magnetoresistance, and the exchange of structural, optical, magnetic, electronic 
characteristics? !^?, The Perovskites (ABX3), A and B are cations and X is an anion. There аге 
many kinds of the perovskite structure crystals, such as the oxides-perovskites (АВОз), flouro- 
perovskites (ABF3), anti-perovskite (A3BX), double-perovskite (А2ВОб) and mnitrides-anti- 
perovskites (A3BN)!*. Furthermore, the perovskite single crystals are considered as proper 
candidates in solar cells, photodetectors, LED, and laser applications because of their facile 
fabrication and excellent optoelectronic characteristics. 


Barium Zirconate, BaZrOs is a surface effect ceramic material that has many mechanical and 
technological applications!é. Also, BaZrO3 has some unique properties such as high dielectric 
constant (^50). In addition, it has high thermal stability and wide bandgap energy (~4 eV). These 
diverse properties have placed this perovskite ceramic as a useful material in applications like 
microwave devices, and high-temperature proton conductor for electric ceramics ?. King-Smith 
and Vanderbilt calculated lattice constants, band structures, zone-center phonon frequencies, and 
elastic constants using the FLAPW method, with the exchange local density approximation 
(РА). Barium Hafnate BaHfOs, which is a compound that has high-temperature melting 
(2893K) and unique physical properties, has been studied in both experimental and theoretical 
terms including thermal expansion coefficient, heat capacity, thermal conductivity, hardness, by 
using the X-ray diffraction pattern, and the lattice parameter of 0.4171 nm at room temperature". 
Also, Vali studied the properties such as the electronic band structure, phonon dispersion and 
dielectric constants of BaHfOs using DFT. He found that BaHfOs has an indirect bandgap 
between (В-Г) points about 3.533 eV". 


Thin films of BaHfi.TiOs were deposited using pulsed liquid injection MOCVD technique on 
Si (100) substrate. The roughness average value (Ra) of the film was increasing with increase x- 
value. Also, the BaHfo.7sTio2s03 thin film showed  well-behavior capacitance-voltage 
characteristic and its k-value is around 1419. Evarestov investigated the structural stability for 
both BaZrOs and ВаН Оз perovskites using hybrid functional DFT with phonon frequency 
calculations. The Mulliken and Born atomic charges results have emphasized the stability of 
cubic and covalent bonds in both crystals!?. 


Photoluminescence and photocatalytic activity of bismuth doped BaZrO3 were investigated 
experimentally using powder X-ray diffraction (XRD), UV-VIS diffuse reflectance 
spectroscopy, photoluminescence (PL) spectroscopy, and photocatalytic activity. The results 
suggested that BaZrO3:Bi could be a good candidate to be used as a visible-light activated 
photocatalyst under excitation wavelengths less than 800 nm”. Parida et al synthesized BaZrO3 
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ceramic powder by the solid-state reaction method. Their measurement indicated that BaZrOs 
could be used for small dielectric resonator antenna (DRA) technology because of its low fewer 
microwave dielectrics materials?'. Khan and Qureshi observed that Tantalum doped BaZrO3 at 
Ta concentration х = 0.04 (BaZro.96Tao.0403) gives the highest rate of hydrogen evolution 
(180umol/h) along with the surface area”. 


Zhang et al conclude that the spontaneous polarization, born effective charge and dielectric of 
BaZrOs is significantly sensitive to epitaxial strain and the transition from paraelectric to 
ferroelectric happened in both compressive and tensile strain. BaZrOs could be useful for tunable 
dielectric materials”. Chen et al synthesis BaZrOs hollow nanostructure with oxygen vacancies 
through one step solvothermal method using ethylenediamine (EDA) as solvent for 
photocatalytic hydrogen evolution from pure water. The H» production rate is about 15 times 
than that of the commercial BaZrO3 (BZO-C)™*. The BaZrO3 hollow nanocrystal dopant by Бе"? 
ions were prepared through facile solvothermal method for photocatalytic Н> from pure water 
under visible light. The existence of Fe (Ш) ion plays a key role in broadening the optical 
absorption range and restricting the recombination of photo-induced carriers under visible light 
irradiation" 220, 


Although these two perovskites are studied intensively, the literature review indicates that 
there is a paucity of thermoelectric data and calculations. The main goal of this research study is 
to review and investigate in details the ferroelectric and thermoelectric properties of the cubic 
Barium Zirconate and Barium Hafnate perovskites to design new potential applications in 
ferroelectric and thermoelectric technologies. 


This contribution is organized as follows: an introduction with the previous overview, the 
computational and theoretical framework details, the results and discussion with the comparison 
with previously available data, and finally the conclusions are drawn. 


2. Computational Details 


We have performed total energy calculations for the oxides-perovskite BaTMOs, an ideal 
cubic structure that contains only one, composed of four atoms in its unit cell. The cations Ba are 
at the cube corners, Та (0, 0, 0), while the cations TM are at the center of the cube, 1b (1, 12, 12) 
and oxygen atoms at the center of the cube faces forming a regular octahedron, Зс (0, У», 2) sites 
of Wyckoff coordinates, as shown in Figurel. The space group for these two perovskites is Pm- 
3m (#221). We have used the interstitial region to express that plane waves are used 
approximation for the calculation of exchange-correlation energy functional, namely; generalized 
gradient approximation (GGA) of Perdew, Burke and Ernzerhof (PBE) parameterization’, using 
modified Becke Johnson exchange potential approximation as parameterized by Tran and Blaha 
(TB-mBJ) to calculate the electronic band structure and density of states??. The TB-mBJ method 
was known that it provides more accurate energy band gap values than PBE-GGA method with 
cheap calculations cost. The first-principles calculations were performed based on (DFT) to 
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describe the ground state??? used the method of (FP-LAPW)*! where this is performed using the 
computer code in Wien2k code**. Also, we chose the interstitial area of the plane wave with a 
cutoff Kmax X Кит = 8, where Кит is the smallest muffin-tin radius in the unit cell, also Ктах 
is the largest K vector magnitude of the interstitial plane-wave expansion of the basis set. The 
chosen number of k-points is (16x16x16) in codes, and which corresponds to the total number of 
k-points of 4096. 


For the optical properties, we calculated the complex dielectric function e(w) = €1(@) + i€2(@), 
where €1(@) is the real part and £»(o) is the imaginary part. The real part ¢1(@) of dielectric 
function e(w) gives from Ше Kramer-Kronig relationship: 





2 оо WE2(0) г 
в1(0) = 1+ = Р 1, at d [1] 
where the Р is the principal value of the integral. The imaginary part ¢2(@) was calculated from 
the momentum matrix elements between the occupied and unoccupied wave functions within the 
selection rules 
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е (0) = (È) Yew f d? k(ckIp"|vk)(vk|pf|ck) x 5( eck - evk — w) [2] 

where p is Ше momentum matrix element between states of bands a and В with crystal 
momentum К. ск and vy are the crystal wave functions corresponding to the conduction, and 
valance bands with crystal wave vector k. There are two contributions to ¢(@), namely; inter- 
band and intra-band transitions. The contribution inter-band is important only for metals, where 
be inter-band between valence-to-conduction absorption spectra. Whereas intra-band between 
valence-to-valence and conduction-to-conduction absorption. Which is describes the complete 
response of a material to the applied electromagnetic radiations field?. Then, the main optical 
parameters: the refractive index п'(о), the extinction coefficient k“(w), the absorption 
coefficient a" (а), optical conductivity oa? (а), electron loss energy function LË (w), and the 


normal incidence reflectivity RË (w) in the crystal are given by***>. 
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[5] 


Веоб(о) = Ime [6] 
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19 (w) = є (ш) + eU? (ш) - e (о) [7] 
В (во) = тті [8] 


The transport properties are obtained from Ше analysis of band structure results based оп 
Boltzmann theory method under the relaxation time approximation as implemented in BoltzTraP 
code?9?7, To explore the properties of thermoelectric such as Seebeck coefficient (or coefficient 
of thermopower) (S), electrical conductivity (c), specific heat (Ce), electronic thermal 
conductivity (K^, and power factor rate (PF) a dense k-mesh of 25000 in the full Brillouin zone 
was used to guarantee good results. 


The Seebeck coefficient (S) is related to carrier concentration via the Mott formula as 
follows: 


ps T?kgT E dn(E) Le] 
— зе ln ав и dE 


[9] 


where е is the electron charge, kg is Boltzmann constant, h is the Planck constant, и is carrier 
mobility and п is the carrier concentration. 


Electrical conductivity (о) is related to carrier concentration as follows??: 


с = [10] 





where т is the electron mass, n is the carrier concentration апа т is the relaxation time. 


Electronic specific heat (Ce) is related to temperature and chemical potential as follows": 


Co(T; W = fnG)G з i0 [PEP de [11] 


3. Results and Discussion 


3.1 Structural Properties 


Goldschmidt tolerance factor (t) for perovskite is given Бу": 


— (ratrx) 
С V2(rg+rx) [12] 


Where ryis the ionic radius of atom N. The t values of existence perovskite compounds range is 
(0.75 « t € 1). For ideal cubic perovskite structure, t has unity value. The tolerance factor t for 
BaZrO3 and ВаНЮз are 1.0 and 1.016, respectively. 
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We presented Ше relationship between Ше energy changes as a function of Ше unit cell 
volume. Using the Birch-Murnaghan equation of state to determine the different important 
structural and ground-state properties?. The equilibrium lattice constants (а) and ground-state 
energies (Eo) are calculated, as estimated from optimization at zero pressure. Besides the volume, 
the lattice constants increase because of the increase in the atomic radius from Zr to Hf. In the 
present work, we were fitting the energy-volume values to a third-order Birch-Murnaghan 
equation of state to calculate lattice parameter a and bulk modulus Bo of Ba(Zr, Hf)O3 
compounds, where the bulk modulus increases as the volume decreasing with moving of cation 
down in the group in the periodic table. The calculated ground-state parameters for both BaZrOs 
and BaHfOs are tabulated in Table 1 along with previous computational and experimental 
studies. For both perovskites, our results are in reasonable agreement with the available 
experimental data and previous calculations. 


3.2 Electronic properties 


To study the electronic properties of BaZrO3 and BaHfO3, DOS’s and band structures along 
the symmetry direction of the first Brillouin zone, of these compounds both PBE-GGA, and TB- 
mBJ methods for exchange-correlation potential are used. The band structures for BaZrOs and 
BaHfOs are presented in figure 2. The Fermi level Er is chosen to locate at 0 eV. In both 
methods, the BaZrOs shows semiconductor characteristics with indirect bandgap (Eg) (А > Г), 
but in the TB-mBJ method, the (Ez) is wider. The valence band maximum (VBM) occurs along 
the R-point symmetry line, while the conduction band minimum (CBM) occurs along the I -point 
symmetry line. The bandgap approaches 3.39 eV by using the PBE-GGA method whereas it 
approaches 4.424 eV by using the TB-mBJ method. As well, BaHfOs shows the same behavior 
with indirect semiconductor bandgap (R > Г), where the bandgap is 3.69 eV by using the PBE- 
GGA method and 5.246 eV when using the TB-mBJ method. Table 2, the calculated bandgap 
(Eg) of both perovskites BaZrOs and BaHfO3 are compared with previous theoretical and 
experimental data. For both compounds, the TB-mBJ values approach enough to experimental 
results. 


The total density of state (TDOS) and partial density of state (PDOS) site projected for both 
compounds are calculated using both methods PBE-GGA and TB-mBJ. The TDOS and PDOS 
for BaZrOs are shown in figure 3a. It is clear that the lower valence band consist of Ba-p state 
around -10.1 to -11 eV, and the upper valence band extends down to approximately -4 eV and is 
composed mainly of O-p states with a small contribution from Zr-d states, while the lower 
conduction band is mainly a mixture of Ba-d and Zr-d states, with a small admixture from O-p. 
In the conduction bands, the Zr-d states are dominant. Figure 3b displays the TDOS and PDOS 
for BaHfOs perovskite using both PBE-GGA and TB-mBJ methods. The contribution to the 
lowest valence band at - 10 eV is due to Ba-p state, whereas the upper valence band (around -4.2 
to 0 eV) is from O-p states hybridized with some Hf-p and НЕд electrons. These O-p and Hf-d 
hybridizations suggest covalent bonding contributions. The bands above the Fermi level are the 
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conduction band of Ше compound; Ше lower part of Ше band (around 3.8 to 6 eV) is dominated 
by the O-p. Also, it can be seen from the plot that the top of conduction bands consist mainly Ba- 
d and Hf-d states hybridized with some O-p. The conduction bands are composed mostly of 
antibonding between Hf-d and O-p as well as antibonding between Ba-d and O-p, maximum 
peak at 7.5 eV is due to the unoccupied 5d states of Hf (5d-eg) antibonding with O-2p mixed 
with Ba-5d and O-2p. 


3.3 Optical Properties 


In this subsection, the optical characteristics of BaHfO3 and BaZrOs are presented. Electronic 
excitations will take place within the crystal if it is exposed to electromagnetic radiation. As a 
result, the optical features of the crystal will be observed. The nature of the response to an 
incident electromagnetic wave at a particular frequency determines if this material could be used 
in optoelectronic devices. The different optical characteristics of the crystals are mainly 
calculated by using the frequency-dependent complex dielectric tensor ғ%(а)). Since BaZrOs and 
BaHfOs perovskites has cubic structure symmetry which means one tensor component is enough 
to do optical calculations. Thus, the optical characteristics of the materials n(w) , k(w), R(o) , 
0 (o) , a(w), and L(w) could be derived from their =; (о) апа €2(w) values as illustrated in the 
computational details section". 


The imaginary part £5(«) is determined from the electronic structure through the joint density 
of states and the momentum matrix elements between the occupied and the unoccupied wave 
functions within the selection rules. In figure 4, the optical characteristics for both BaZrOs and 
BaHfO3 compounds are plotted. In panel a) the spectrum of €,(w) shows the three peaks 
structure for BaZrO3 and BaHfO3. The (w) has deeply impact on the absorption spectrum of 
the matter. A large value of =, (о) refers to a large absorption spectrum include of the whole of 
all different transitions that will occur between the valance and the conduction bands. It is 
starting to response to the incident radiation at critical points which are the threshold energy at 
2.87 eV апа 4.26 eV for BaZrOs and BaHfOs respectively. These points are C,-C. splitting and 
give a threshold for indirect optical transitions between the highest valance band and the lowest 
conduction band. This is known as the fundamental absorption edge. The main peak in the 
spectrum is situated at 7.12 eV and 7.52 eV for BaZrOs, and BaHfOs, respectively. These peaks 
are dominated by transitions from the O-p band and the small contribution of transition metals d 
state just below Fermi energy to Ba-p state of the conduction band. In panel b) the spectrum of 
€1(w) for BaZrO3 апа BaHfO3 compounds show that the static dielectric constant ғ1(0) are 
equal 3.69 and 3.51, whereas the main peaks appear at energy 5.24 eV and 5.59 eV, respectively. 
It is important to consider both sı and e2 together since they affect each other, meaning the shape 
of є» cause corresponding changes in the shape of є; and vice-versa. This is known as the 
Kramers-Kronig relation between the real (с1) and imaginary (62) parts of the dielectric function. 
In short, the dielectric function describes what an electric field such as an oscillating light wave 
does to the material. In panel c) the absorption coefficients for BaZrO3 and BaHfO3 compounds 
show the same behavior as є. («). The absorption edges start from 3.28 апа 4.29 eV for BaZrO3 
and BaHfOs respectively which are corresponding to optical bandgaps. This produces from the 
transition of the O-p valence band electrons to the empty d conduction bands in Zr or Hf. In 
ВаНГО: perovskite, the main peak in dielectric absorption spectra occurs at 7.81 eV whereas the 
two other peaks at 11.3 and 11.78 eV whereas in BaZrOs the main peak occurs at 7.3] eV the 
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other two peaks occur at 8.8 and 11.0 eV. In both compounds, Ше Ba does not contribute to Ше 
low-energy absorption spectra; because its electron valence bands are not in the transition energy 
domain. In panel d) the loss function for BaZrOs and BaHfO3 compounds are plotted and their 
peaks occur at 7.6 eV and 8.4 eV, respectively. In panel e) the reflectivity for both BaZrO3 and 
BaHfO3 compounds are plotted and their static reflectivity R(0) is approximately equal 10%. At 
low energies (1.е., 0—7.8 eV for BaHfO3 and 0-7.2 eV for BaZrO;), the reflectivity spectrum was 
observed to extend as high as 33% and 30% for BaHfO3 and BaZrOs;, respectively. Our 
calculations show a strong reflectivity at high energies increase up to 75%, and 72% for BaHfOs, 
and BaZrOs, respectively. In panel f) the refractive index plots show that the static refractive 
index n(0) are equal 1.92 and 1.87 for BaZrOs and BaHfO;, respectively. The refractive index 
reaches its maximum value 2.7 at 7.5 eV for BaHfOs, whereas for BaZrO3 its 2.5 occurs at 5.0 
eV. Finally, the optical conductivity o(w) for both BaZrOs and BaHfO3 perovskites are 
displayed in panel g). There is quite similarity in curves of both perovskites. For BaZrOs the first 
peak appears at 6.8 eV, the second peak appears at 9.0 eV whereas the third peak at 10.5 eV. On 
the other hand, at BaHfOs there аге two mean peaks at 7.4 eV апа 11.5 eV. The maximum value 
for BaZrOs is equal 8900 at 10.5 eV while for BaHfOs it is equal 9100 at 11.5 eV. 


3.4 Thermoelectric (TE) Properties 


The efficient thermoelectric characteristics for any materials determines by their figure of 
merit ZT = S*oT/k, the increasing of S and с increase figure of merit while increasing the total 
thermal conductivity к decrease figure of merit. As it is known, a suitable material for 
thermoelectric applications should have figure of merit equal or greater than one which give 
energy transformation efficiency greater than 2590 at optimal operating temperature". 


Figure 5a shows the variation of Seebeck coefficient S with respect to temperature for both 
compounds BaZrO3 and BaHfOs. The two curves have approximately the same behavior where 
they change nonlinearly by increasing from zero to reach the maximum value 206.86 uV/K at 
T=270K for BaZrO3 and 288.13 un V/K at T=180K for ВаН Оз, after that it is decreasing slightly 
in ВаН Юз but in BaZrOs it is still increasing. The anomalies behavior of Ше Seebeck coefficient 
at the low temperature 15 due to the strong electron-phonon coupling, where the phonon drag 
effect becomes manifest. The peak in the Seebeck curve around 100 K in the figure is coming 
from the phonon drag effect. It is roughly consistent with the estimated value of 05/5, where др 
being the Debye temperature of both compounds [05 (BaZrOs) = 544 К and 05 (BaHfO3) = 465 
К]. At room temperature (RT), Seebeck coefficients for BaZrO3 and BaHfOs are 208 апа 261 
uV/K, respectively, which are close to conventional thermoelectric material values. The positive 
value of S in the given temperature range defines the majority charge carriers as holes i.e. p-type 
semiconductors. High S implies the large effective mass due to the existence of energy bands 
near the Fermi level in this material. In figure 5b, the electric conductivity o as a function of 
temperature when the constant relaxation time Tọ = 0.8 х 107145 is used, it is shown that the 
electric conductivity is decreasing non-linearly at the low temperature to reach its global 
minimum at T=121K in BaZrOs case and T=90 іп BaHfOs, then it starts increasing linearly with 
temperature increases for both BaZrO3 and BaHfOs. The electric conductivity at RT for BaZrOs 
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is about 0.554 x 1055/т and it about 0.289 x 1055 / т for BaHfO3. Figure 5c displays Ше 
electronic thermal conductivity x? as a function of temperature when the constant relaxation 
time то = 0.8 х 107115 is used. The behavior of the electronic thermal conductivity curves are 
nonlinear increasing for both compounds. The electronic thermal conductivity of BaHfOs at RT 
is equal to 0.154 Wm'!K'!, which is about half the BaZrOs value of 0.309 Ут КЛ. Our results 
contradict the previous experimental values 10.4? and 5.24 Wm''!K"!, respectively. This 
contradiction could be because our calculations are based on Boltzmann transport theory which 
unfortunately calculates only the electronic part (x?) of total thermal conductivity (x = к? + Kr), 
whereas, the experimental results measured the total thermal conductivity (x), where the lattice 
part (ку) contributed as the dominant part in both compounds. The ratio between electronic 


0 
thermal conductivity and electronic conductivity obey Wiedemann-Franz law c — LT), where 


L is the Lorentz number (in order of ^10 92WOQK-? W) for both compounds. The temperature 
dependence of the thermal conductivities of BaZrO3 and BaHfOs follow an inverse temperature 
law (1/T) at low temperatures, showing typical phonon conduction characteristics ^9. Figure 5c 
shows how are the power factor (PF) curves changed when the temperature increase for both 
compounds. They reach their maximum values 0.213 uW/mk? at T=180K for ВаН Оз and 
0.0814 uW/mk? at Т-240К for BaZrOs. Also, at RT the PF's equal 0.128 un W/mk? for BaHfOs 
and 0.0778 uW/mk? for BaZrOs. Although the PF values for both compounds are small 
compared with the well-known thermoelectric materials, these results support that BaHfOs is 
better than BaZrOs if a suitable element dopant in specific concentrations is used and/or through 
band engineering or switching to nanostructures. It is possible to fabricate perovskite at the 
nanoscale as the lattice thermal conductivity decreases. This will improve the merit figure at 
room temperature; consequently, the thermoelectric efficiency increases. 


To determine in details thermoelectric (TE) behavior for both BaHfOsand BaZrO3, such as 
Seebeck coefficient (5), electrical conductivity (с), power factor (PF), electronic specific heat 
(Се), and electronic thermal conductivity (X^) have been expressed as a function of chemical 
energy (и); the calculated characteristics were plotted for considered temperatures 
range: 300-800 K. In figure 6, S for both Ва/гОз and BaHfO3 compounds as a function of 
chemical energy и at different temperature are plotted. For both compounds, the curve increases 
linearly from zero value to reach the maximum value and then decreases sharply to negative 
global minimum values; after that it increases linearly again to reach zero value. 


The chemical potential dependencies of electric conductivity rate (с/т) for both BaZrOs and 
BaHfO3 compounds at different temperatures are depicted in figure 7 panels a) and b). In the two 
panels, the electric conductivity goes to zero in energy range 0.52 — 0.74 Ry for BaZrOs and 0.54 
— 0.8 Ry for BaHfOs for the wide range of temperatures from 200K to 800K. In figure 8 the 
power factor with respect to chemical energy at different temperatures is plotted; a) BaZrO3 and 
b) ВаНЮ: perovskites. The two plots (a and b) display two main peaks one for n-type and the 
other for p -type. The p-type of BaZrOs exhibits the PF peak 11.7 x10!! W/mK?s at 0.525 Ry, 
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while Ше n-type 6.0 x10!! W/mK?s at 0.83 Ry. For BaHfOs, the p-type and n-type PF peaks 
10.65 x10!! W/mK?s and 6.75 x10!! W/mK?s are at 0.5 Ry and 0.775 Ry, respectively. 
Therefore, these materials may be used as p-type thermoelectric materials. Also, the PF values of 
the two compounds around Ев are negligible РЕ(и = Ер) = беу. This means that their 
thermoelectric properties are not good without doping. Figure 9 shows the electronic specific 
heat capacity as a function of chemical energy at different temperature. The calculated specific 
heat capacity of BaZrOs in this study is well agreement with the data in the literature". The C, 
values for ВаН ЮЗ are slightly greater than that of BaZrOs because of the higher molar mass of 
Hf compare to Zr. Figure 10 shows the variation of the density of state around the Fermi level 
with respect of temperature. For the two alloys BaZrOs and BaHfOs the band gap still appear for 
a wide range of temperature. 


4. Conclusions 


Employing the all-electron full-potential LAPW method based on DFT, with using PBE-GGA 
functional of exchange-correlation energy; and semi-local exchange potential TB-mBJ to study 
physical properties of both BaZrOs and BaHfO3 perovskites. Also, semi-classical Boltzmann 
transport theory has been used to calculate the thermoelectric properties. A summary of results 
are: 


(i) The calculated lattice constants at zero pressure of these compounds are in a reasonable 
agreement with the previous theoretical and experimental values. 


(ii) The electronic bandstructure calculations showed that BaZrO3 and BaHfOs are insulators 


materials with an indirect bandgap (А - Г). The fundamental bandgap of this compounds 
increases when mBJ-GGA approach is used instead of PBE-GGA. 


(ii1) The joint DOS and the Kramers-kronig transformation are used to calculate and analyze the 
optical characteristics: dielectric constant e(w), the absorption coefficienta(w), Ше 
reflectivity R(w), the refractive index n(w) and electron loss energy L(w). The oxygen p states 
and Transition metals Zr- and Hf-d states play a major role in optical transitions. 


(iv) Finally, the calculated transport characteristics showed that the two compounds are not 
suitable choice for thermoelectric applications but if they are dopant with some impurities their 
thermoelectric properties could be improved. 
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Table caption 


Табіе!: Lattice constant а(А), bulk modulus B(GPa), for Ше perovskites BaZrO3 and BaHfOs. 








Compound a(A) B(GPa) Method 
4.234 149.8 This work 
4.1895, 4.19348, 4.199749. 4. 19254650, 1034550 155.98, Exp 
Ва7гОз | 4.22698, 4.2063°!, 4.18122!, 4.181572, 4.19326 1271! 
4.1547,4.20733, 4.14817, 4.257“, 4.193"), 175.37,157, 174.7", | Theory 
4.162329, 4.192, 4.19810 14111, 1722, 16810 
4.212 156.0 This work 
4.1714645 4.1756, 4.1727", 4.16685, 4.17963, Exp 
PE 4.1711, 4.192, 4.19310, 4.310%, 4.24999, 1443. 1752, 1712. Тпеогу 
4.17018, 4.2867, 4.24861, 4.16082, 4.0649? 123.6”, 186", 194% 
146.72!" 145.501?, 
1798 








да, 128 Obtained from Ше Voigt-Reuss-Hill (VRH) approximation. 


95. 125 Obtained from fitting the third-order Birch-Murnaghan equation 
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Table 2: The calculated bandgap (Eg) using PBE-GGA and TB-mBJ methods with comparison 
with previous computational and experimental results 








BaZrO; BaHfO; 
Eg (eV) method Ref. E, (eV) method Ref. 
4.42 TB-mBJ This work 5.25 TB-mBJ This work 
3.39 PBE-GGA This work 3.69 PBE-GGA This work 
5.3 Exp ol 6 Exp (predicted) 18 
3.9 LDA di 2.99 LDA Е 
3.156 ССА 5 3.94 SX-LDA 12 
5.4 LCAO 10 3.17 GGA и, 
3.55 Ехр. Е 5.3 SX-LDA 93 
4.96 Exp. 25 3.9 ССА-РВЕ 52 
5.9 TB-mBJ Е 
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Figure Captions: 


Figure 1: Crystal structures of ВаТМО: (TM: blue, Ba: green and О: red). 

Figure 2: Band structure of cubic perovskite: а) Ва/гОз PBE-GGA, b) BaZrO3 TB-mBJ, с) 
BaHfO; PBE-GGA, and b) BaHfO; TB-mBJ methods. 

Figure 3: Calculated total and Partial-projected DOS: a) Ва/гОз using PBE-GGA compared 
with TB-mBJ and b) BaHfO3 using PBE-GGA compared with TB-mBJ method. 

Figure 4: Calculated optical properties for the BaHfO3 and BaZrO3 compounds, а) €2(w), b) 
#1 (о)), c) a(w), d) L(w), е) К(о), f) n(o), and в) с(о). 

Figure 5: Transition coefficients as а function of temperature at ground state chemical potential 
for the BaHfO3 and BaZrOs compounds: а) Seebeck coefficient (S), b) electric conductivity rate 
(с/т), c) electronic thermal conductivity (ко), d) The power factor (РЕ). 

Figure 6: Seebeck coefficient S as a function of chemical potential at different temperature for 
a) BaHfO3 and b) BaZrOs. 

Figure 7: (a) and (b) are the electrical conductivity rate (с/т) as a function of chemical potential 
at different temperature for both BaHfO3 and Ва/гОЗз; (c) and (d) are the electronic thermal 
conductivity as a function of chemical potential at different temperature for both compounds 
BaHfO; and BaZrOs, respectively. 

Figure 8: The power factor (PF = 576/т) аз a function of chemical potential at different 
temperature range for а) BaHfOs and b) BaZrOs. 

Figure 9: The electronic specific heat as a function of chemical potential at different temperature 
for a) BaHfO3 and b) BaZrO3. 

Figurel0: The density of state at Fermi level DOS [Er] as a function of chemical potential at 
different temperature for a) BaHfOs and b) BaZrOs. 
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